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Role of shuffles and atomic disorder in Ni-Mn-Ga 
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We report results of ab-imtio calculations of the ferromagnetic Heusler alloy Ni-Mn-Ga. Particu- 
lar emphasis is placed on the stability of the low temperature tetragonal structure with c/a = 0.94. 
This structure cannot be derived from the parent L2i structure by a simple homogeneous strain 
associated with the soft elastic constant C'. In order to stabilise the tetragonal phase, one has 
to take into account shuffles of atoms, which form a wave-like pattern of atomic displacements 
with a well defined period (modulation). While the modulation is related to the soft acoustic 
[110]-TA2 phonon mode observed in Ni2MnGa, we obtain additional atomic shuffles, which are 
related to acoustic-optical coupling of the phonons in Ni2MnGa. In addition, we have simulated 
an off-stoichiometric systems, in which 25 % of Mn atoms are replaced by Ni. The energy of this 
structure also exhibits a local minimum at c/a — 0.94. This allows us to conclude that both shuffles 
and atomic disorder stabilize the c/a = 0.94 structure. In both cases the stability seems to be 
associated with a dip in the minority-spin density of states (DOS) at the Fermi level, being related 
to the formation of hybrid states of Ni-d and Ga-p minority-spin orbitals. 

Keywords: Ni-Mn-Ga alloys, martensitic phase, modulated structure 
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Ni-Mn-Ga alloys (close to stoichiometric Ni2MnGa) 
are known to exhibit unique magneto-elastic properties. 
They are ferromagnetic at room temperature (Tq ~ 
380 K) , and undergo (including the precursor) a two-step 
martensitic transformation for Tm < Tc (Tm ~ 200 K) 
[1]. In the martensitic state of Ni-Mn-Ga the structure 
consists of differently oriented martensitic domains (twin 
variants), which are also magnetic domains. This makes 
the martensitic structure of Ni-Mn-Ga sensible to an ex- 
ternal magnetic field, which can induce a redistribution 
of the martensitic domains in the sample. Those domains 
with easy magnetic axes are along the field will gain in 
energy on cost of the domains with less favourable ori- 
entation of the magnetisation. Alignment of twin vari- 
ants by the motion of twin boundaries can result in large 
macroscopic strain up to 6 % [2, 3]. 

This effect is used in the magnetic-shape-memory 
(MSM) technology [4]. The MSM technology is based 
on the magnetic field induced redistribution of marten- 
sitic domains in the sample. From a technological point 
of view, Ni-Mn-Ga is more promising than other ma- 
terials being presently in commercial use, for example, 
the well-known material Tb-Dy-Fe (Terfenol-D) which 
exhibits magnetostrictive strains of about 0.1 %. De- 
sign of new efficient MSM magneto-mechanical actuator 
devices is in progress [5] . 

In this work we discuss the stability of different struc- 
tures inside a single martensitic variant of Ni-Mn-Ga. 
From experimental studies we know that these alloys 
can form at least three different phases in the marten- 
sitic state [6]. Depending mostly on composition, crys- 
tals can be found in the modulated tetragonal structure 
with c/a = 0.94 known as 5M or lOM, the orthorhom- 
bic modulated structure 7M, or the tetragonal structure 



with c/a ~ 1.2. Sometimes, the modulated structures 
are also denoted by 5R and 7R [7]. 

The central question concerns the nature of the mod- 
ulation in the 5M and 7M structures, which has often 
not often taken into account in theoretical investigations 
[8-10]. However, the modulation plays a basic role in 
Ni2MnGa, as shown by recent first-principles calculations 
[11]. In the literature, it is argued that the reason for the 
modulation arises from specific nesting properties of the 
Fermi surface in Ni2MnGa, which in turn, causes soften- 
ing of the acoustic [110]-TA2 phonon mode [12]. 

In the previous study [11] of the stability of the 5M 
structure, we have drawn attention to specific changes in 
the DOS, namely a dip, which develops in the minority- 
spin density of states (DOS) at the Fermi level, being 
related to the Fermi surface geometry [13]. Analysis of 
the partial contributions to the DOS shows that this dip 
is formed by Ni-dJ, and Ga-p| electrons. Hence, there 
is a strong evidence that the formation of alike hybrid 
states might contribute to the stability of the modulated 
structures. In addition to the shuffles leading to the mod- 
ulation, atomic disorder and non-coUinearity of magnetic 
moments were proposed as possible factors being also re- 
sponsible for the structural stability of the modulated 
phases in Ni-Mn-Ga. 

In this work, we discuss properties of the modulated 
5M structure with new features which have been obtained 
in recent calculations. In addition, we consider an off- 
stoichiometric composition of Ni-Mn-Ga with excess Ni 
replacing Mn atoms. For the total energy calculation 
we use the Vienna Ab-initio Simulation Package (VASP) 
[14,15] and the implemented projector-augmented wave 
formalism (PAW) [16]. The electronic exchange and cor- 
relations are treated by the generalised gradient approx- 
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FIG. 1: Modulated martensitic structure of Ni2MnGa (5M) 
shown schematically: (a) projection of the 5M structure on 
the (001) plane (top view) and (b) projection of the 5M struc- 
ture on the (110) plane (side view). Filled, open and gray dots 
show the positions of Mn, Ni and Ga, respectively. In addition 
to the modulation, which moves the atoms along the [110] di- 
rection only, there is a tiny shuffling of the atoms of the order 
of 0.001 A (~ 1 % of the modulation amplitude), which moves 
the atoms in [001] and [110] directions. The extremely small 
displacements of atoms have been enlarged in the figure for 
clearness. 



imation (GGA). The 3d electrons of Ga have been m- 
cluded as valence states. 

In order to simulate the 5M structure, we used a super- 
cell, which resembles the five-layered structure obtained 
by V. V. Martynov and V. V. Kokorin in experiment 
[6] . The supercell consists of five tetragonal unit cells of 
Ni2MnGa, similar to those used in Ref. [9]. This allows 
us to incorporate the full period of the 5M modulation 
in the supercell [11]. The basal plane is spanned by [110] 
and [001] of the L2i structure with the modulation along 
[110]. Altogether we use ten atomic planes perpendicular 
to [110] in order to form the supercell. The modulation is 
generated by displacing these atomic planes along [110] 
direction. By this construction, two full five-layered pe- 
riods fit into the supercell. The initial magnitudes of the 
atomic displacements were chosen according to the 5M 
structure of Ref. [11], for the Ni-plane equal to 0.324 A 
and for the Ga-Mn-plane equal to 0.292 A. 

This supercell has an orthorhombic symmetry with lat- 
tice parameters a ~ 4.17 A, 6 « 20.73 A, c ~ 5.633 A, 
yielding a tetragonality ratio of c/a w 0.955, i.e, we take 
exactly the structure which was obtained theoretically in 
Ref. [11]. In the calculations we allowed for relaxation 
of the structure by changes of the volume, cell shape 
and atomic positions. The plane-wave cutoff energy was 
equal to 241.6 eV and the k-points were generated using 
the Monkhorst-Pack method with a grid of 10 x 2 x 8 
points in the full Brillouin zone for the long orthorhombic 
supercell. Actually, all parameters of the calculations are 
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FIG. 2: Off-stoichiometric supercell of Ni2+a;Mni_a;Ga for 
X = 0.25 created by replacing Mn by an Ni atom in one of the 
16 positions of the L2i structure. 



the same as used in our previous calculations [11], except 
for the k-points mesh. 

The structure obtained after relaxation is shown in Fig. 
la and b. As in the model for the 51VI structure in Ref. 
[6,11], the modulation has form of a static wave with po- 
larization (the direction in which the atoms move) along 
[lIO], which is perpendicular to the propagation direction 
[110]. All atoms move with the same phase that agrees 
with the acoustic mode [12,17]. However, a careful anal- 
ysis shows that the shuffling of the atoms in Ni2MnGa 
consists of two different contributions, being a superpo- 
sition of them. The first one is the modulation which 
is wave-like. While the second one is different in the 
sense that it is not a wave. We call these additional shuf- 
fles tetrahedral distortions^ similar to those considered by 
Harrison for tetrahedrally coordinated solids having the 
symmetry of diamond or zincblende [18]. According to 
Harrison, radial and angular distortions of tetrahedral 
structures stand for this motion of the atoms. In case 
of Ni2MnGa each Ni atom is surrounded by two tetra- 
hedrals formed by Ga and Mn, respectively. Note that 
amplitudes of the tetrahedral distortions are of the order 
of 0.001 A which is ^1% of the modulation amplitude. 

We would like to draw attention to some features re- 
lated to the tetrahedral distortions. In the (001) plane. 
Fig. la, the Ga and Ni atoms move closer to each other, 
while the Mn atoms move towards the free space left 
by the Ni atoms. Also, the plane (110) in Fig. 1(b) 
shows couples of Ni atoms moving closer to their nearest 
Ga atom, while the nearest to them Mn atom is pushed 
out. Thus, each atom of Ga tends to couple with two 
Ni atoms, which will contribute to the hybridization of 
the Ga-p and Ni-d orbitals as discussed above. We stress 
that the modulation and the tetrahedral distortions must 
be considered separately from each other. The first one 
is related to the soft acoustic [110J-TA2 phonon mode, 
while the second one is connected with the coupling of 
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FIG. 3: The relative change of the total energy of Ni-Mn- 
Ga as a function of cja for the three different cases: cu- 
bic Ni2MnGa structure, modulated 5M structure and off- 
stoichiometric case defined by Ni8+iMn4_iGa4. 



acoustic phonon modes [110]-TAi and LA with optical 
modes derived mostly from the Ni atoms, as it was dis- 
cussed in Ref. [17]. 

As a matter of fact, most experimental studies con- 
centrate on ofF-stoichiometric compositions of Ni-Mn-Ga, 
while theoretical studies are mostly done for for the stoi- 
chiometric structure Ni2MnGa. Hence, theoretical inves- 
tigations for off-stoichiometric compositions are urgently 
needed. The results obtained so far allow to discuss im- 
portance of the Ni-Ga interaction in Ni-Mn-Ga. show- 
ing that the stability of the modulated shuffles benefits 
from this coupling (and the locally broken symmetry). 
The modulation facilitates the formation of hybrid Ga- 
•p states with other atoms, which is less favourable in 
a perfectly symmetric environment, where Ga "cannot 
choose" neighbours to couple (the shuffles of atoms re- 
move this degeneracy allowing for a lower energy of the 
system). 

The question now is whether the modulated shuffles are 
the only way to remove the perfect order. This could also 
be done by doping the structure by local defects. Planing 
to check this we have simulated off-stoichiometric Ni-Mn- 
Ga structures. In order to simulate an off-stoichiometric 
composition, a cubic supercell of 16 atoms is used in the 
calculations. The corresponding structure with a com- 
position formula Ni84_iMn4_iGa4 is shown in Fig. 2. 
An extra atom of Ni is placed into a Mn position (cen- 
tral atom in Fig. 2). The supercell has been optimized, 
first, with respect to the atomic positions and the vol- 
ume. The cja = 1 ratio was kept constant at this stage 
of the calculations. After the relaxation, the optimal 
volume became smaller (192.51 A^) compared to the sto- 



ichiometric structure (195.79 A^) by about 1.6%. The 
atomic positions have changed as well, but only for the 
Ni atoms. The positions of Mn, Ga and the extra-Ni 
remained unchanged, which is not physical, but is due 
to the choice of the supercell and the periodic boundary 
conditions used for the calculations. One would need to 
simulate a bigger supercell for a more sophisticatd treat- 
ment of the defects in Ni-Mn-Ga. In Fig. 2 the extra-Ni 
is shown in the center connected by bonds to its eight 
nearest Ni neighbours. The relaxation has led to consid- 
erable displacements of all eight Ni atoms towards the 
extra-Ni along the bonds, while the rest of the structure 
remains unchanged. The magnetic moments remained 
ferromagnetically aligned. Although, the cxtra-Ni has 
got a smaller magnetic moment (0.24 /is) compared to 
the regular Ni atoms (w 0.35 ^b)- 

This relaxed cubic structure is now used to study the 
impact of a tetragonal deformation. Fig. 3 shows the de- 
pendence of the total energy on the c/a ratio for the three 
cases: the perfect stoichiometric Ni2MnGa structure, the 
5M structure, and the off-stoichiometric structure defined 
by Ni8+iMn4_iGa4. It turns out that the total energy 
curve for the case of the off-stoichiometry has a minimum 
at c/a = 0.94, exactly the value as for the 5M structure, 
and in the experimental investigations. This shows that 
both modulation and disorder are responsible for the sta- 
ble tetragonal structure with c/a = 0.94. 

Fig. 4 presents the total electronic density of states 
for the off-stoichiometric case (structure shown in Fig. 
2). The DOS is shown for: the cubic structure with 
c/a = 1 and the tetragonal structure with c/a = 0.94. 
Important is here the dip in the minority-spin density of 
states right at the Fermi level. The dip is present for 
both cases, but stronger developed for the more stable 
tetragonal structure. The analysis of the partial DOS 
shows again that the two peaks, responsible for the dip, 
originate from the Ni-d and Ga-p orbitals (the effect is 
similar to the results obtained for the 5M structure in 
Ref. [11]). 

In summary, we have presented computational re- 
sults for the modulated 5M structure of Ni2MnGa. We 
find that shuffles in Ni2MnGa involve two different ef- 
fects, which contribute to the pattern of the atomic dis- 
placements. The first and the stronger contribution, 
called modulation, arises from the soft acoustic [110]-TA2 
phonon mode, which is observed in the phonon spectrum 
of Ni2MnGa. The second effect, referred in our work as 
tetrahedral distortions, yields small (as compared to the 
modulation) atomic displacements and results from the 
coupling between the acoustic [110]-TAi and LA with 
corresponding low-energy optical modes of Ni. The cal- 
culations for the off-stoichiometric Ni8+iMn4_iGa4 su- 
percell yield a local minimum of the total energy at 
c/a = 0.94, which is the "natural" tetragonality ratio for 
the Ni-Mn-Ga at low temperatures. The same tetrag- 
onality ratio can be obtained in the calculations when 
taking modulation (5M structure) into account. Simi- 
lar to the case of the 5M structure, the stability of the 
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FIG. 4: Electronic density of states of the off-stoichiometric 
alloy Ni2+i:Mni_a;Ga with x = 0.25 for the two cases: unsta- 
ble c/a = 1 and metastable c/a — 0.94. (a) Total DOS for 
the cell with 16 atoms, (b) Enlargement of the total DOS 
which shows closer the dip in the minority spin density of 
states. Analysis of the partial density of states shows that 
the states which form the dip are derived from the Ni-d and 
Ga-p orbitals. 



c/a = 0.94 ratio is related to a dip in the minority-spin 
electron density of states which develops from hybridiz- 
ing Ni-d and Ga-p states right at the Fermi level, stabi- 
lizing the tetragonal variant with c/a = 0.94. In other 
words, the tetragonal structure becomes stable due to 
the covalent interaction of Ni an Ga atoms via the p-d 
hybrid electronic orbitals, the importance of which for 
the Hcuslcr alloys have been duscusscd by Kiiblcr et al. 
[19]. But we emphasize the role of the local symme- 
try loss which facilitates the formation of the p-d hybrid 
states. In future work we will perform calculations using 
the GGA-I-U method. This will split apart the bonding 
and anti-bonding states around the Fermi level, which 
will further deepen the dip discussed above. We expect 
that the resulting total energy will correspond then to 
the 5M structure as the ground state. 
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